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Abstract. The scattering of an incident plane wave on two Aharonov-Bohm vortices 
with opposite fluxes is considered in detail. The presence of the vortices imposes non- 
trivial boundary conditions for the partial waves on a cut joining the two vortices. 
These conditions result in an infinite system of equations for scattering amplitudes 
between incoming and outgoing partial waves, which can be solved numerically. The 
main focus of the paper is the analytic determination of the scattering amplitude in 
two limits, the small flux limit and the limit of small vortex separation. In the latter 
limit the dominant contribution comes from the 5-wave amplitude. Calculating it, 
however, still requires solving an infinite system of equations, which is achieved by the 
Riemann-Hilbert method. The results agree well with the numerical calculations. 
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1. Introduction 

The Aharonov-Bohm problem certainly leads to one of the most fascinating quantum 
mechanical effects which has been experimentally tested [2]. 

Clearly the original model as stated in the seminal 1959 paper jl] is by itself 
somehow an abstraction. Leaving aside the zero flux size limit issue and the related 
question of the behaviour of wave functions at the location of the vortex (see, e.g., [3]) 
there is not such a physical system where the flux of a magnetic field can pierce a plane 
at a certain location without returning at some other place on the plane so that the 
total flux is zero. In the Aharonov-Bohm paper this issue is solved by assuming that the 
return flux is spread far away on the plane (in fact at infinity) so that its direct effect 
on the problem at hand can be neglected. This simplifying assumption is crucial in 
rendering the model exactly solvable, and most of the research effort has been devoted 
to this situation (see, e.g., |5] and references therein). 

We would like to address here a more realistic situation where the piercing flux 
returns at a finite distance. More precisely, we consider two vortices of opposite strength 
piercing the plane at distance 2R one from the other. Doing so, we move from the 
solvable one-vortex standard Aharonov-Bohm problem to a two-vortex problem, which 
is more difficult to address, with no explicit solution at disposal. 

A few many-vortex Aharonov-Bohm problems have already been considered in the 
literature, in particular the scattering on vortices arranged on a lattice [6] or vortices 
whose locations are random [7j. 

The problem of scattering of an incident plane wave on several vortices has mostly 
been discussed by two methods. In [8]-[in] a formal expansion of the scattering 
amplitude in a series of Mathieu functions has been constructed and in [II]-[T3] a 
formal diagrammatic-like series has been proposed. So far these approaches have had a 
limited success and did not produce explicit formulas except when the distance between 
the vortices goes to infinity [T3] . 

Here we consider the scattering problem on two vortices in the singular gauge with 
non-trivial boundary conditions along a cut connecting the vortices. To address this 
problem, we use an approach already put forward in [H] for studying the low-energy 
spectrum of a charged particle in an harmonic well coupled to two Aharonov-Bohm fluxes 
of different strength. In the plane a discontinuity has to materialize on any branch cut 
joining the two vortices due to the non-trivial Aharonov-Bohm phases accumulated by 
the charged particle moving around one or the other vortex. The positions of the vortices 
are fixed at R and —R on the x-axis, so that the branch cut can be chosen to be the half- 
circle of radius R centered at the origin in the lower half-plane (of course the observables, 
here the modulus squared of the scattering amplitude, should not depend on the choice 
of a branch cut). As here the total fiux is zero, no phase is accumulated at infinity 
by the particle encircling both vortices. The scattering amplitude is, by definition, a 
series of scattering amplitudes between incoming and outgoing Hankel partial waves. 
The non-trivial boundary conditions on the cut lead to an infinite system of equations 
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Figure 1. The branch cut in the plane and the boundary conditions for two Aharonov- 
Bohm vortices with fluxes a and —a. 



for the expansion coefficients, which are solved numerically. The main part of the paper 
is devoted to an analytic investigation of two limiting cases. The ffist is the small flux 
limit (more precisely, the small a limit, where a is the dimensionless ratio of the flux 
to the quantum of flux) and the second is the small R limit (in fact the small kR limit, 
where k is the scattering momentum). In both cases the analytic expressions obtained 
are in good agreement with the numerical calculations. 

The plan of the paper is the following. We start in Section |2] by developing the 
formalism for the scattering of an incident plane wave making an angle 9 with the x- 
axis. We illustrate in Section [3] the physics at hand with several numerical simulations 
in various cases of interest. Next, in Section H] we consider the small a limit and 
find an exact lowest-order expression for the scattering amplitude. In Section [5] we 
consider the small kR limit, which allows for some simplifications in the contribution 
of the relevant partial wave scattering amplitudes. Then in Section |6] we transform the 
small kR equations into a form which lets us find an exact expression for the scattering 
amplitude by the Riemann-Hilbert method. In Section [7] we show how the same solution 
can be obtained and simplified by applying the Riemann-Hilbert method directly to the 
problem. In Section |H] we outline the main results. A few details of the calculations are 
given in the Appendix. 

2. General formalism 

We consider two Aharonov-Bohm vortices on a plane at points R and — i?, i.e., with 
polar coordinates (-R, 0) and (-R, vr) with fluxes a and —a respectively (see figure [I]). 
We work in the singular gauge, where the electromagnetic potential is removed but the 
wave function is defined with a cut connecting the two vortices such that on the opposite 
branches of the cut the wave function and its normal derivative obey the relations 




(1) 
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and 

^¥'\RA) = f{<P)§^^^'-\R,<P). (2) 

Here and below superscripts (+) and (— ) indicate functions outside and inside the circle, 
respectively, and /(0) describes the necessary jump along the cut. In the geometry of 
figure [1] it can be written in various equivalent forms 

^^^^~\ 1, 0^/ "^1 0, 0^/-' \e-™ 0^/ 

where / = [vr . . . 27r] denotes the angular interval of the cut. 

Everywhere except the cut the wave function obeys a free Schrodinger equation 

(A + A;2) ^ = 0. (4) 

Inside the cut circle the wave function can be expanded into series of regular elementary 
solutions of (jl]). We find it convenient to use the Bessel function basis 

oo 

¥-\kr,^)= CmJ\m\{kry'^t (5) 

m=— oo 

Let us denote the wave function outside the circle by \I'(+)(/cr, 0). To implement the 
boundary conditions ([1]) and ([2]), one multiplies both equations by e"""^ and integrates 
over from to 27r for all integer m from — oo to oo. One concludes that 

J{m\ (^) j^^ e-™^vl>(+) (x, 0)/(0)d0 - J|„| (x) j^^ e-''"'^^ (vI/(+) (x, 0)) /(0)d0 = (6) 
where x = kR. 

The usual way to exploit these conditions (cf. [H]) is to expand \E'("''^(r, 0) into a 
series of Hankel functions provided that the incident wave is fixed: 

oo 

^;+)(A;r,0) = J|i|(A;r)e"<^+ ^ t„(/)i/|5(A;r)e^"'^. (7) 

n=— oo 

From ([6]) it follows that the coefficients t„ (/) obey an infinite system of equations 

oo 

("^Iml-f^H ^ ^H^lnl'^ Amntnii) = -{J\m\J\l\ " J\m\J\l\)Aml , (8) 

n=— oo 

where the AmnS are the Fourier coefficients of /(0) in (j3]) 

A^. = ^ r H4>y'--'m = e"" f J!r-^.-™ I" = "I . (9) 
/vr /n — — ^ m ^ n 

Here and below, when the argument of Bessel functions is not specified, it is meant to 
be equal to kR = x. It is convenient to rescale the tn{l)'s as 

tn{l) = (10) 



Scattering on two Aharonov-Bohm vortices with opposite fluxes 

so that the ?/„(/) 's obey the system of equations 

T' H^^^'\ f V V 




. , A^^y^iS) = -\^-^\A^,. (11) 

Jm Hn' I \Jm Jl / 

From figure [1] it is clear that the scattering problem is symmetric with respect to 
a refiection in the ?/-axis. Namely, if \l/(r, 0) is a solution then \l/(r, vr — 0) is also a 
solution. This symmetry is a consequence of the invariance of the function under 
the transformation 

/(7r-0) = /(0). (12) 

It manifests itself in a symmetry of the matrix elements A^n 

^—m,-n — ( 1) ^mn (l^) 

and of the expansion coefficients ?/„(/) 

y_„(-/) = (-l)"+'i/„(/). (14) 

The matrix Amn has two more symmetries. A trivial one corresponds to an integer shift 
a 1 + a, which leaves Amn and y„(/) invariant: 

?/„(/; a) = ?/„(/; 1 + a). (15) 

This invariance refiects the well known fact that the Aharonov-Bohm effect depends 
only on the fractional part of a. 

The second symmetry is related with the transformation a ^ —a. From ([9]) one 
concludes that 

A^^{-a) = e-2™A_„_„(a) = (-l)"^-"e-2™A„„(a) (16) 

which leads to the symmetry for the coefficients ?/n(0 

y„(/; -a) = y_„(-/; a) = (-l)"+'2/„(/; a). (17) 

The knowledge of coefficients t„(/) (or y„(/)) determines all other quantities. Of 
particular interest is the amplitude of scattering of a plane wave on the two vortices. 
To determine it, one has to find a solution which at large distances from the vortices is 
the superposition of an incoming plane wave in the direction 6 and a circular symmetric 
outgoing wave 



T--s>oo V T^kr 

Expanding the incoming plane wave into a series of Bessel functions (see, e.g., [T5] ) 

oo 

gi/crcos(e-0) _ il'lj|^|(A;r)e"('^-^) (19) 

Z=— oo 

the scattering wave function can be expressed through \l/['''^(A;r, 0) defined in ([7]) as 

oo 

^(+)(fcr, 0) = XI '^^^^^'f\kr, 0)e-"^ = e^'^'^^^^^"'^) + m^"'^\kr, 0) . (20) 



l=—oo 
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the reflected field given by a series of Hankel functions of the first kind 

oo 

vl/M)(A;r,0)= i'"tn(O^H • (21) 

n,i=— oo 

From the asymptotic behaviour of the Hankel functions (see, e.g., jl5] ) 



iJ«(A;r) ^ J^e'^'^™/'-"/^) (22) 

r-5-oo V Tlkr 

one gets that the scattering amplitude F{6, (p) in f|T8l) is 

oo 

F{e,(f))= J2 i'''''"'tn(Oe''"^"'''' • (23) 

n,/=— oo 

The symmetries f[T^ and f[T7|) lead to the symmetries of the scattering amplitude 

F{n-e,n~(l)) = F{e,(f)) (24) 

and 

F{9, 0; -a) = F{9, (p;! - a) = F{n + 9,7r + 0; a) . (25) 
These relations make it possible to restrict the value of flux to < a < |. 

3. Numerical results 

We proceed to calculate the scattering amplitude numerically, by solving the systems 
(fTTj) for yn{l), for different values of /, and substituting the result into (fTO|) and then 
(1231) . From physical considerations it follows that contributions of partial waves with 
|/| -C kR have to be small. In the calculations below we restrict ourselves to |/| < lOkR 
and check that higher waves do not change the results noticeably. 

Truncation in n is more complicated. In order to find the coefficients themselves 
with good precision, one has to retain much higher values of n in the systems (fTTj) . 
\n\ < N with N ~ 100. Even then, the precision attained at technically feasible values 
of N is not sufficient, and one has to extrapolate the results to — oo. 

For a given finite N, we truncate the infinite sum in ( fTTl) to n = — + 1, . . . , A^, 
so that the system, for a given /, consists of 2A^ equations for as many variables. 

It turns out that convergence in is improved by employing the technique used 
in [H], discretizing the boundary conditions. The issue at hand is that by choosing 2A^ 
coefficients t„(/) in ([7]) and, respectively, 2A^ coefficients in ([5]), one can, in general, 
satisfy the boundary conditions ([I])-(l2l) exactly for no more than 2A^ discrete values 
of 0, rather than for any 0. A convenient recourse is to enforce those conditions at 
0fc = (^~|);§^) ^ = I)---; 2A^ (these 2A^ points are distributed uniformly on the circle 
and avoid as much as possible the locations of the fluxes, where the wave function is 
singular). Then, integration over in ([6]) has to be replaced with summation over 0^. It 
is easy to see that the only change this entails is a modified expression for Amn, namely 

cosTTo; (m = n) 

^mn=e'™<( siriTra l-C-l)"-'" / „\ • (26) 
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For — > oo, one recovers OH]), whereas for a given finite N , the resulting amphtude is 
closer to the infinite- one than the one obtained without this modification. 

Convergence in of the resulting scattering amplitude is slow, due to the usual 
problem inherent in Aharonov-Bohm or anyon numerics: We are trying to represent a 
wave function which has a fractional power behaviour at the positions of the vortices 
(cf. f lll6p ). as an expansion in terms of regular wave functions. The corresponding 
coefficients fall off slowly. We have found, however, that for all values of kR, 6, and a 
considered, and for any G [— vr, tt], the A^ dependence of the amplitude fits an empirical 
formula 

F(e.^;N)^F(e,4>) + ^ + j^+'^ + j^ (27) 

well, even for rather small A^. (We have used the values A^ = 40 . . . 100 for the fitting.) 

In figures [2}|3] we present the results of numerical calculations for the amplitude of 
the scattering on two Aharonov-Bohm vortices for fluxes a = 1/4, 1/3, and 1/2 with 
kR = 1 and with incident angles 6' = 0, 7r/4, 7r/2. For clarity, in these figures we count 
the refiection angle from the direction of the incident wave, i.e., we plot F{6,6 + 0) 
instead of F{9, 0). 

A notable feature of the results is the absence of symmetry with respect to the 
refiection — )■ —0. This is in contrast to the standard Aharonov-Bohm scattering, 
where, despite the fact that chiral symmetry is broken by the presence of the vortex, 
the amplitude remains symmetric under such refiection. One could have expected the 
same here for ^ = 0, when the geometry of the system is just as symmetric with respect 
to the oncoming wave as it is with a single vortex. However, in general this turns out not 
to be the case because for the invariance the spacial refiection has to be accompanied by 
interchanging of the vortices. For = only the case a = 1/2 is symmetric with respect 
to inversion — )■ —0 (cf. and ( 125|) ). In general, with positive a the refiected wave 
is more likely to deviate to the right (negative 0). Only for 6 = tt/2 is the amplitude 
symmetric, as anticipated (cf. 



4. Small a limit 



To find the leading term of the scattering amplitude in the limit a — )■ 0, it is convenient 
first to transform the boundary conditions to another form which is of independent 
interest. Multiplying both sides of (E]) by il™le'"^(^'+'^), summing over m, and using again 
the identity f lT9|) . one concludes that has to obey 



2tt 



ixcos(^?'-0)_^ / ^(+)(X, 0)) - ¥+\x, 0)^ f ^-■^-cosie'-^)\ 



d 



/(0)d0 = 



(2J 



This equation can also be obtained by a direct application of the Green identity inside 
the circle. 




Figure 2. Modulus squared of the scattering amplitude for kR = 1 and the incident 
angle 9 = 0. Left: a = 1/4, center: a = 1/3, right: a = 1/2. 
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Figure 3. The same as in figure [2] but for the incident angle 9 = 7r/4. 



For the pure scattering setup one has to take \l/*^"'')(A;r, 0) as in (120|) . From ( 128|) it 
follows that the reflected field obeys the equation 



2-IT 



j{e,e') 



/(0)d0 



where 



27r 



d 



d 



doc 



dx 



-ix cos{9' ~ij>)\ \xcos{0 



(29) 

/(0)d0. (30) 
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Figure 4. The same as in figure [2]but for but for tlie incident angle 6 = 7r/2. 



The expression in the square brackets is a total derivative: 

ixcos(6»'-(/))_^ (^ixcos{e-<f>)\ ^ /g-ixcos{e'- 

dx ' dx \ 

e + e' 



e 
2i 



-^\x cos{6—(f)) 



, 9-e'\ d 



—2ix sin 



sm 



6-6' 



so the integration in fl29]) for given in Q easily follows as 



4 . 



J(0, 0') = — e^™ sin vra cot sin ( 2x sin 



6-6' 6 + 6' 
sin 



X 2 V 2 2 

Using the second expression for in ([3]), fl29|) can be rewritten as 



^6,6') = -2ie™sin7ra 
d 



2tt 



dx ' 



d0 



4 e-^' / 0-^' 6 + 6' 

H e sm na cot sm 2x sm sm 

X 2 V 2 2 



where 



(31) 



(32) 



(33) 



d0 . (34) 



y Q Ox Ox V / 

It is well known (and easily checked either by the Green representation of the reflected 
field or by the direct substitution of f l2T|) in flM|) ) that ^{6,6') is proportional to the 
scattering amplitude for the reflected field 



^{6,6') = -F{6,6') 



(35) 




Figure 5. Scattering amplitudes on two vortices witli a ~ .01 and kR = 1 for 
different incident angle Q versus the reflection angle Left: 6* = 0, center: B — 7r/4, 
right: = 7r/2. Solid lines are direct numerical solutions of pT|) . Dashed lines are 
calculated from the analytic formula (|37)). 



If one assumes that at small a (hence, small sinvra) the wave function ^!'^^^^\r^(j)) can 
be expanded into a series in a as 

^M)(r, (j)) = a^i(r, 0) + a^^2{r, (j)) + . . . (36) 

(the zero order term is absent since when a is an integer there is no Aharonov-Bohm 
scattering) and takes into account that the right-hand side of fl33l) is proportional to 
sinvra, then the scattering amplitude in the leading order can be approximated by 

F(0, ^') ^ — i sin vra cot — ^ — sin (2x sin — ^ — sin — ^ — j . (37) 

The existence of a term proportional to A^~" in the extrapolation formula ( |27|) renders 
numerical calculations as very small a with of the order of a few hundred uncertain. 
Fortunately, the analytical formula ( 1371) for the dominant contribution at small a (1371) 
is regular and does not contain singularities at the vortex positions. This means that 
the coefficient in front of A^~° is, at least, of the second order in a. Therefore, at very 
small a we drop that term from the fit (!27|) . 

The results obtained in this manner are presented in figure [5] for different incident 
angles and compared with fl37|) for small a. For clarity, we divide the scattering 
amplitude by sin vra and, like in the previous Section, count the reflection angle from the 
direction of the incident wave. The agreement is quite good and the differences between 
both curves are of the order of the terms which have been neglected (i.e., sinvro;). 

5. Small kR limit 



Another interesting limit corresponds to the vortices close to each other, so that kR <^ 1. 
In this case it is convenient to use the rescaled variables (fTTl) and to expand all quantities 
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at small x. One gets (see, e.g., [15]) 



Jln{x) \m\ H^a^'{x) -\n\+6noP 



Jra{x) X ' hI^\x) 



X 



(3J 



with 



^ 7r + 2i(ln(x/2)+7) ^^^^ 



(7 is the Euler constant). 

In this limit equations (ITT]) take the form 



^ (|m| + \n\ - 5nop)Amnyn{l) = - \l\)Ami (40) 

n=— 00 

SO that the dependence on a; = kR is entirely contained in p. On the other hand, from 
( |T0|) it follows that when x — )• „(/). Therefore at small x the dominant 

contribution to the scattering amplitude comes from the S'-wave amplitude (/ = 0), 
more precisely from yo{0), as it is usual for the scattering on small-size objects (see, 
e.g., [IE]). Nevertheless, to determine yo{0), one still has to solve an infinite system of 
equations PUI) with / = 0: 
00 

^ (|m| + \n\ - 5nop)^mnyn{0) = -\m\Amo ■ (41) 

n=— oo 

The scattering amplitude (j23ll is dominated at small x by yo{0) 

F{e, <P) to(0) = ^ -i^pyo(O) (42) 

where we used that, at small x, Hq^\x) ~ 2i/(7rp) with p given by ( l39l) (see, e.g., [T5]). 
From ffT4l) it follows that 



y_„(0) = (-l)"i/„(0) (43) 

and one can transform (1411) to a system of equations involving the ?/„(0)'s with non- 
negative n only. When m = one obtains 



oo 



2 sinvra / n 

cosTra pi/o V(l - (-l)")l/n(0) = . (44) 

71 ^ — ' 

n=l 

For positive m, straightforward calculations demonstrate that due to symmetry ( H3|) 

oo 

{m + \n\)Aranyn{0) 

n=—oo, n^O 

oo . oo 

= 2m V A„„y„(0) + ^^(1 + (-1)-) V(l - (-l)")y„(0) . (45) 

IT ^-^ 

n=l n=l 
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Using (Hlj) and adding the n = terms, one gets that (HTl) with positive m is equivalent 
to 

sin vrafl — (—V 



^ ] ^m,nl/n(0) 



n=l 



27rm 
/ sinvrafl — (— 1) 



:i+i/o(o)) 

cosTrafI + (— 1) 



+ 



27rm2 ' 4m 

All the dependence on x = kR in this equation is in p only. 

It is convenient to redefine the variables ?/n(0) with > 1 as 

VniO) = ^x„(l+i/o(0)) 

and to set 

Xo = 1 

so that fH6l) takes the form 



n=0 



where 



and 



Similarly 



^ ^ A-mri'^n ^frr 

2pyo(0) 



2' 



1 + 1/0(0) 

sin 7ra(l — (— 1) 
27rm^ 



cos7ra(l + (-1)™) 



4m 



becomes 

00 

-^cot7ra = ^x„(l-(-l)") . 

n=l 

The system ( H9ll is linear, so its solution has the form 

where a„ and 6„ depend on a but not on p and not on yo{0). From 
that ^ has to be determined by 

00 00 
cot na = Y, a„(l - (-1)") + ^ - (-1)" • 



n=l 



n=l 



Solving for ^ one finds 

P2/o(0) 1 



1 + 2/0(0) /3(a) 
where P{a) is real. Therefore 

1 



Z/o(0) 



l + /3(a)p 

where the dependence on kR and a has separated. 



mio). (46) 



(47) 
(48) 

(49) 

(50) 
(51) 

(52) 

(53) 

one concludes 

(54) 

(55) 



(56) 
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From fl42p the scattering amplitude in the small kR limit is 

fc^o ~7r + 2i(ln(A;i?/2)+7 + /3(a)) ' ^^^"^ 

It has the standard form of a short-range scattering amplitude (e.g., the scattering on 
a small disk of radius R with Dirichlet boundary conditions corresponding to P{a) = 
(see, e.g., [IE])). 

A simple way of determining P{a) numerically is to fix p to a given value po (say, 
Po = 0.3) and to approximate the infinite system (HOl) by a truncated one 

N 

^ (|m| + \n\ - Snop)Amnyn{0) = -\m\Amo (58) 

n=-N 

with N finite (m is also taken from — to N). fl58l) is solved numerically and yields 
?/o(0) for a given as a function of a, from which /3(a) follows using fl56l) . The true 
/3(a) is then obtained by taking the limit A^ — )■ oo. As the resulting function should 
have power singularities at the vortex locations (cf. (11161) below), the convergence with 
increasing A^ is slow, all the more so when a is close to 0. We found that at large A^ a 
good fit of the numerical data involves the same terms as in ([2] 



CLl 0,2 O,^ CL4, /rr\\ 

For illustration, in figure E] we present on the left numerical calculations for a = .2 and 
p = .3 for A^ from 150 to 550. The best fit to the data corresponds to (159|) with 

ao - -2.2688, ai ~ -.4504, 02 ~ -133.8, 03 ^ 449.3, ~ -552.0 .(60) 

On the right of figure El the difference between the data and the best fit is presented. 
It is small (~ 10^^), structureless, and can be attributed to random numerical errors, a 
strong argument for the validity of the fit. Of course, this type of extrapolation is not 
stable, depends on the form of the extrapolating curve, and the results should be taken 
with some care. 

In figure [7] we present numerical estimates of /3{a) for a few values of a calculated 
by the above method together with the exact expression fl89|) which will be derived in 
the next Sections. 

In figure |8] the scattering amplitudes calculated numerically directly from the main 
system of equations (ITT]) for kR = 0.01 and kR = 0.1 are plotted for a = 1/4, 1/3, 1/2 
and the incident angle 6 = 0. (For other incident angles the pictures are similar.) The 
dashed lines in this figure represent the small kR prediction f lSTj) . where we used the 
following values of /3(a) computed from the exact formula 



/3Q^ =-21n2, /3Q^ =ln2-^ln3, /3Q^ =-ln2. (61) 

The agreement between numerics and the asymptotic expression is good, with differences 
between them being of the order of kR, as expected. 
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Figure 6. On the left: red circles are the values of ?/o(0) calculated numerically from 
(|4T]) with a — .2 and p = .3; the dashed line represents the limiting value of ?/o(0) (i-e. 
oq in the fit ([59])). On the right: the difference between yo(0) computed numerically 
at finite N as on the left figure and the fit ([59]) with fitting coefficients (|60| . The solid 
red line serves only to join the points. 





Figure 8. Black solid lines: modulus squared of the scattering amplitude with the 
incident angle 9 — and different fluxes. From bottom to top: a = 1/4, a — 1/3, and 
a = 1/2. Blue dashed lines are the predictions ([57]) for the corresponding values of a. 
Left: kR = 0.01, right: kR = 0.1. 



6. Exact solution for small kR 

The system of equations ( H9l) permits a useful interpretation. Let us define the 
generating functions of x„ and /„ 

oo 

<^^^\z) = J2^nz\ (62) 

n=0 

and 

oo 

F{z) = J2fnz'' (63) 

n=l 



which rewrites, using ( jSTj) . as 

^, , cosvra ,^ \ , /. \x sinvra,^. , , ^. , ,, 
i^(^) = —— (ln(l - ^) + ln(l + z)) - —— (Li2(2) - Li2(-2)) (64) 
4 zvr 

(Li2(-2) = Xl^i -2"/"^^ is Euler's dilogarithm). The normalization f l48|) leads to $*^+)(0) = 

1. 

In principle, the coefficients of an expansion of a function holomorphic 

outside the circle (see Section [7]). It would be more appropriate to use instead of 
z = re"^ the variable 1/z = e^'^/r. To simplify the notations, we use the variable z, but 
retain the superscript (+) in the function fl62|) to stress that it is related to a function 
well behaved outside the circle. Equation ( H9l) can then be interpreted as 

<l>W(e''^)/(0) = ^^'\e''^) + ^F{e''f') (65) 
where ^^~\z) has an expansion in non-positive powers of z 

$(-)(.) = 5:^. (66) 

n=0 
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Equation f l65p is the standard equation for the non-homogeneous Riemann-Hilbert 
problem (see, e.g., jl^)' ^^"^ solution can be obtained by the usual methods [T7] . 
First one needs to find the solution of the homogeneous equation 

TW(e''^)/(0) = T(-)(e'*) . (67) 
which can be done by taking its logarithm 

lnT(+) -lnT(-) = -ln/(0) . (68) 
An elementary solution of the last equation is the Cauchy integral [17] 



hiT{z) = -- 



2n 



ln/(0) 



■d0 



(69) 



2tt Jq 1 - ze-"t 

Here for symmetry reasons we use the third definition of /(</>) in ([3]) and again drop the 
factor e^'^", hke in fH9|l : 



m = 

It follows that 

tW(^) = 

One easily checks that 



< < TT 

TT < d) <2n 



1 + 


)■ 










cot- 









, T^-\z) 



1 + 1/2 
1 - l/z 



< < TT 
TT < d) <2'K 



(70) 



(71) 



(72) 



so that dSID is fulfilled. 

Using (1671) . the full non- homogeneous equation ( l65l) can now be rewritten 

T(-)(e''^ 



as 



${+)(ei'/ 



and therefore 



$(+)( 



T(+)(ei<^) 
) <l>(-)(e'' 



$(-)(e''^) +ei^(e''^) 
) ^ , i^(e''^) 



T(+)(e''^) T(-)(ei^) ^T(-)(ei^) ' 
A solution of ( FMl) is again given by a Cauchy integral 

r(2) 27r io T(-)(ei</')[l-2e-i'^] ^ ■ 

Note that since -^(2;) and T^^\z) are regular respectively inside and outside 
circle, the integral fl75|) cannot be calculated as residues. In the next Section 
that it can be simplified by using a different approach (cf. (I113P ). 
A general solution of (!75|) has the form |17] 

^{z)=iT{z)J{z)+T{z)P{z) 

where P{z) is a polynomial and 



tan" (0/2)- 



F(e'<^ 



+ e 



-ma/2 



ze~ 



cot"(0/2) 



F(-e'<^ 
1 + ze- 



(73) 



(74) 



(75) 

the unit 
we show 



(76) 



(77) 
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Under — vr — (I77|) has the symmetry 

Jiz) = J{z) . (78) 
In particular, for z real J{z) is real and can be written as 



1 

J{z) = - Re 

vr In 



tan" (0/2) + e-'™/2 cot" (0/2) ^ 



, . d0 . (79) 

Using $*^+)(0) = 1 and the fact that ^^^\z) starts at infinity with a constant we infer 
from (|76l) that 

^{z) = iT{z)[J{z) - Jm+T{z) . (80) 

Clearly (IHOj) is a solution of (!65|) and consequently of fH9l) for all ^'s. Eq. (152!) has yet 
to be fulfilled. Formally it can be rewritten as 

-^cot7ra = lim[$(l- 5) -$(-!)] (81) 

2 (5-j-O 

where ^{z) is the same function as in fl52]) which is given by f lSU]) . ($(1) in general 
diverges which explains the necessity of the indicated limit.) 
According to (lA.lip and (lA.20p of Appendix, one has 

$(-!) = -^cot(7ra)e (82) 



and 



"^(1 - = ( i ) Wi^) - -m) + 1] + ^ cot(vra)e + 0{5'-n . (83) 



Comparing these behaviours with flHTj) . we conclude that ^ has to be chosen in such a 
way that the singularity of fISU]) at 2; = 1 cancels 

e(J(l) - J(0)) + 1 = 0. (84) 

It means that the required solution of f H9|) is 



^/ ^ ^ J(l) - J(Z) 

^^")=^(")j(TW(o) 

where T{z) and J(z) are defined in (!7T|) and (177|) respectively. Expanding $(2;) into 
power of z gives the solution of (H9l) . 

One can now calculate all quantities of interest, in particular, P{a) in fl55|) 

/3(«) = 2(J(1) - J(0)) , (86) 

so that using (1791) one finally gets 



P(a) = - Re e^""/'tan"(0/2)— ^-e-^""/^cot°(0/2) .\ ■ 

Note that under a — t- 1 — a, F{z) — )■ —F{—z), the integrand in the last expression 
remains symmetric and therefore 

(3{1 - a) = (3{a) , (88) 

as it follows from 



7) 
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Using (11131) . one can prove that when < a < 1 (1^71) is equal to 

/3(a) = \^ ) + \^ (1^) + 7 + In 4 + — ^ . (89) 
2 \2J 2 \ 2 / 2 sm vra 

Here the function "^{z) is the logarithmic derivative of the gamma function and 

7 = —^1/(1) is the Euler constant. /3(a) calculated numerically from this expression 

is presented in figure [3 



7. Riemann-Hilbert approach 

In the previous Sections we started from a system of equations fHT!) valid for small 
X = kR and derived an exact solution using the Riemann-Hilbert method which 
matches two holomorphic functions regular inside and outside a given contour. Here we 
demonstrate that the same solution naturally arises from general considerations without 
explicit transformations of boundary equations. 

When kR 0, one can ignore in the vicinity of the vortices the k"^ term in the 
Schrodinger equation so that the wave function obeys a Laplace equation. It means 
that it is a sum of a function of the variable z = re^"^ and a function of the variable 
z = re~^'^. Due to the symmetry ( !T4l) . the S'-wave function should necessarily have the 
form 

^(r,0) = $(2) + <l(-^) . (90) 
Inside the circle the function ^{z) has to be regular and can be expanded into a series 



m z 



00 

¥'-^^^^\z) = J2<^m{^y. (91) 



.R 

m=0 



Outside the circle the situation is more complicated. The expansion with the 
definition f lTOj) states that the outside function has the form 

n=-oo H'n'{kR) 



When kR — )■ one gets 



^[+) {kr, (f)) = (z) + {-z) + pyo{0) In - (93) 

R 



with p given in f l39p . Here 



n=l V ^ / 



(94) 



The logarithmic term appears because of the presence of H^\kr) whose short-distance 
behaviour requires to add ln(r/i?) to the class of admissible functions. \n{r / R) is neither 
analytic nor anti-analytic but can be written as a sum of functions with the same 
symmetry as in (1901) 



In ( ^) = i In f-i^) + J In f i^ ) . (95) 



RJ 2 \ RJ 2 \ R 
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This argument leads to the conclusion that outside the circle the allowed functions have 
the form 



|,(outside)^^) = <|(+)(^) + ^P2/o(0)ln (-i^ 



(96) 



The function $(+^(2;) is regular outside the circle but the logarithmic function is not. 
In the splitting (pOj) it is implicit that if the function ^{z) has a cut at z = \r\e^'^° with 
fixed it should also have a cut at z = |r|e'*^'^~'^°''. Therefore one has to define ln(— i^;) 
on the unit circle z = e^*^ with two cuts, the first one from 1 to 00 and the second one 
from —00 to —1, 



f < < TT 

3tt 



(97) 



2 TT < < 27r 

Now the Aharonov-Bohm boundary conditions ([T]) are reduced to the condition that the 
inside and outside functions of variable z (and of —z) are related by 

|,(inside)(^gi<^^) = /(0)<|(°"*^''i'=)(e''^) . (98) 

Let us first find two holomorphic functions F^^\z) and F^~\z) regular, respectively, 
outside and inside the circle such that their difference obey 

F(-\z)-F^^\z)=f{<P)g{<P) (99) 

with /(0) given in ([3]) and g{(j)) in (1971) . 

The explicit form of F{z) is given again by the Cauchy integral [17] 

Using for /(0) the expression (170|) one gets 



F{z) 



where 



1 

An 



7r/2 



+ 



ze~ 



2-IT 



ze~ 



cos vra Ji + sin Tra J2 (101) 



and 

Therefore for \z\ < 1 
F^~-\z) 
which coincides with 



ze 



ze~ 



1 — 2;e 



cosvro; 



1 + ze- 



d0 



-[\n{l + z) + \n{l- z)] (102) 



l.[U2{z)-Lh{-z)] . (103) 



sm na 



^ [ln(l + ^) + ln(l - z)] - 2^ 
and for > 1 
F^+\z) =F^-\-l/z) . 
Knowing F{z) allows to rewrite (1981) as 

|.(-)(e'<^) = /(0)<|.(+)(e^*) - F(+)(e'^)p|/o(0) 



[Li2(^)-Li2(-^)] (104) 
(105) 
(106) 
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where 

^(-\z) = |.(i^«idc)(^) ^ F'^-\z)pyo{0) (107) 

is a function regular inside the circle and ^^'^\z) defined in (iQlI) is regular outside the 
circle. 

To transform this equation to the form solved in the previous Section, we note that 
due to the invariance (fT2l) /(tt — 0) = /(0) and fll06p remains valid under — > tt — 

|,(-)(_e-i<A) = j(0)|,(+)(_e-i*) _ F(+)(_e-i^)p2/o(o) . (108) 
According to (HI) ?/_„(0) = (-l)"yn(0), so that 

$(+)(_e-i<^) = l(l+^,(0))$(e'*) (109) 

with $(-2) defined in ( l62l) . Denoting <l>*^^)(— e"^"^) by <l>'^^)(e~^'^) and taking into account 
that 

(_e-i<^) = F(-) (e''^) = F(e''^) (110) 
with F{z) given in flM|) . we conclude that fllOSp takes the form 

$(+)(e'<^)/(0) = <l>(-)(e'*) + ei^(e^*) (111) 

which coincides with fl65|) . 

Imposing the condition that the wave functions remain finite at both vortices will 
fix the value of yo{0) in the same way as above (cf. fl84|) ). Therefore in the end we get 
the same solution as in the previous Section. 

But this approach gives us more than just recalculating the solution by a slightly 
different method. From ( l98|) together with ( l96l) it is clear that the problem under 
consideration is equivalent to 

where the definition fp7|) has been used. 

This equation is also a non-homogeneous Riemann-Hilbert equation of the same 
type as (l65l) with substitution F(e''^) — > f{(j))g{(j)). It means that its solution is the 
same as above but the integral (1771) is replaced by 



1 



.i,,/, tan"(0/2) w2S^n0/2) 



j(z) = — e"'™/" vr^/^ V2Z_^ (113) 

To prove that this expression and fl77|) are the same we notice that according to fl99l) 
and ([lOlD 

F[e"f) - F(-e-'*) = f{<P)g{<P) (114) 

so F(e"'^) — f{(f))g{4>) contains only negative powers of e'''^. When this difference is 
substituted into the integral ( ffSll the contour can be shifted to infinity and, therefore, 
the integrals f lll3p are equal to f l77p . Other equivalent forms of this integral will be 
discussed elsewhere. 
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A few general comments about the obtained solution are in order. Combining all 
terms together, we find that in a vicinity of the vortices but outside the cut circle the 
S'-wave function in the limit kR — )■ takes the form 

¥+\r) = ^(1 + |/o(0)) \^{z) + ^{-z)\ + p2/o(0) In l^] (115) 

where r has coordinates x and y , z = {x + iy) / R, z = {x — iy)/R. 

In this expression, ^{z) is given by (185|) with z — )■ 1/z. The scattering amplitude 
yo{0) from the H'^\kr) contribution in the expansion (I7l) is fixed by ( |56l) and ( !87l) . 
That value of |/o(0) renders the function $(^) finite at both vortex locations. In fact, 
as a consequence of fl82l) and (l83l) . the total wave function flllSp vanishes at the vortex 
positions, as it is the case for the Aharonov-Bohm effect 

^ I 1^— -RT, r close to the right vortex with fiux a /_ , 

\I'(r) ~ < ^ ^ ^ . (116) 

I |r + r* close to the left vortex with fiux — a 

In deriving these limits, we take into account that due to the necessary substitution 
z — )■ 1/z, singularities at z = ±1 in ( l85l) and in flllSp are interchanged. 

In this Section we use only one of the two Aharonov-Bohm boundary conditions 
( l98l) . which matches the function on the cut circle, but do not discuss explicitly the 
second condition (|2]), which matches the derivative. The reason is that the functions we 
consider are either analytic or ant i- analytic. But for any such function the derivatives 
with respect to r and are proportional to each other: 

or r 0(p or r 0(p 

Therefore, when condition (p8|) is satisfied and the function /(0) is piecewise constant 
(like in our case) the derivative with respect to r will obey the same condition as the 
function itself, except for 5-function contributions at the points of discontinuities of /(0). 
One can check that (jUj) is related with the cancellation of these discontinuities. As has 
been demonstrated in the previous Section, that equation leads to the vanishing of the 
wave function at singular points (cf. ( I116P ). so the 5-functions give no contributions and 
the derivatives over r automatically obey the correct boundary condition (jj]). 

In [18] it was argued that a solution of the S'-wave type which remains finite at 
the two vortex locations cannot exist. But a logarithmic term like the one in (IllSp was 
not properly taken into account in that discussion. When ignored (which means setting 
^ = in ( 185|) ). only the singular solution T{z) is obtained. The role of the logarithmic 
term is precisely to remove the singularity at the vortex locations. 



8. Conclusion 



We have considered the scattering problem on two Aharonov-Bohm vortices with 
opposite fiuxes a and —a, separated by a distance 2R. First, we developed a numerical 
method for the construction of scattering amplitude. In its simplest formulation it 
reduces to solving an infinite system of equations f lTTj) . The convergence of finite 
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approximants is slow due to the singular behaviour of the wave function at the vortex 
positions (cf. flll6p ). Then we analytically constructed the dominant contributions to 
the scattering amplitude in two particular cases: small values of a and small values of 
kR. 

In the small a limit, i.e., when sinvro; — ?■ 0, the scattering amplitude takes the form 
F(9,9') — > — isinvra cot sin I 2x sin sin ) (118) 

sin-Tra— >0 2 V 2 2 / 

where 6 and 6' are respectively the incident and reflection angles. 

In the small kR limit, the solution has been obtained by the Riemann-Hilbert 
method in its simplest setting. The dominant contribution is given by the S'-wave 
scattering amplitude, which has the form 

^'^ k~^o ~7r + 2i(ln(A;i?/2) + 7 + /3(a)) ^^^^^ 

where (3{a) is 

/3(a) = i^f-) + +7 + ln4 + — (120) 

' 2 V2/ 2 V 2 / ' 2sin7ra ^ ^ 

with "^{z) being the logarithmic derivative of the gamma function. 

The same method permits also to find the behaviour of the S'-wave function in the 

vicinity of the vortices. Denoting this wave function as 

¥^\r-^) = Mkr) + UkR) -J^^H^)[kry-^ (121) 

n=_^Hn [kR) 

we found that when kR — > 

1 + p[a)p 

where 

_ 2i 

^ ~ TT + 2i(ln(A;i?/2) + 7) ' ^^^^^ 
The generating function of y„(0) = y-n(O) in this limit takes the form 

E^"(0)^" = ^'^(^)-2 • ^^^^^ 



n=0 



Here 



and 



'^W = ^(^) W)-Jiz)] (125) 



1 



tan"(0/2) „i^,/2 Cot"(0/2) 



TT 



j(z) = — e"^™/' ^ZZ^ + e-°/^ U__ d0. (126) 

A discussion of other partial waves in the limit kR — )■ and of related questions 
in the general case of two Aharonov-Bohm vortices with arbitrary fluxes will be given 
elsewhere. The exact solution for the scattering amplitude with kR finite remains an 
open challenge. 
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To calculate the behaviour of the solution ( IHOll at 2; — )■ 1 and z — )■ —1, one has to 
investigate the corresponding limits for the integral ( 179|) . When z = —1 + 6 and 5 — 0, 
the second term in ( frf|) is formally diverging. This divergence appears at small (p and 
can be isolated by splitting the integral ( 1791) into two parts: 



J(-l + (5) 



1 



TT 



tt/2 



Re 



>"/2tan"(0/2)- 



Ffe^' 



+ e-^™/2cot°(0/2) 



l + (-l + 5)e-i 



(-1 + 5)6-^"^ 
d0 



7r/2 



Re 



-iW2 cot" (0/2) T 



i(0-i5) 



Changing the variable in the last integral — )■ t(5 gives 
j(-l + 6)-^W{6) + 0{l) 

(5->0 

where 

1 /•'r/(25) 

= - / Re 
Jo 

From known expressions 

Li2(l) = ^ , Li2(-1 



+ 0{1). (A.l) 



(A.2) 



cot"(^t/2) ; ^ 



i(t - i) 



dt 



TT , , TT 

— , Lis -1 = 

it follows that the function F(— e^"^) has the following expansion 



F(-e" 

Finally, one has in the limit 5 — )■ 

W{6) -- 



1 TT TT 

> - COS TTo; m 2(1) — i— H — sinvra 
4 ^^2^8 



Hi Re 



00 g-i7ra/2 



i7rt°(t 



1 n / e \ ''''n ^ 

- cosvra ln(2ot) + i— H — sinTra 
4 ^ ^ ^ 2^ 8 



dt 



The first integral in this expression can be calculated by contour integration 



t°'{t — i) sin Tia 

and the second one by taking the derivative with respect to a of this result 
""^ dtlnt TT^ i 



+ cot vra 



i(:"(t-i) sinvra 2 
Substituting these expressions into (1A.6P and taking the real part, one obtains 



and 



iy(5) = - 



J{-l + 5) 



(5->0 



TT cot(7ra) 
4 



TT cotfvra) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(Ai 



(A.9) 



(A.IO) 
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TC 



$(-1 + S) = -- cot(7ra)^ + 0{5' 
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(A.ll) 



The behaviour of J{z) close to z = 1 can be obtained by a similar method. The difference 
is that in this case the value J(l) is finite but due to the singularity of T{z) at z = 1 
we also need to know the correction term. One has 



J(l -6) = -Re 



■k/2 



^i7ra/2 ^ a 



tan"(0/2) 



+ e 



-i7rQ/2 .a 



COt"(0/2)- 



1 - (1 - 5)e-^^ 
F(-e'<^) 



1 + (1 - S)e-'^ 
First we split this integral into two parts as follows 



d(f) 



Tr/2 



d0 



. . . d0 + 



7r/2 



(A.12) 



(A.13) 



where 5 ^ e ^ 1. 

In the first integral, the first term dominates at small 0, and 



J(l-<5)— > J(l) + r(<5) 

(5— s-O 



where 



V{S) 



1 



Re 



TT 



^i7ra/2 ^ a 



tan'^{(f)/2)F{e"^ 



iS) 10 



d(j) 



(A.14) 



(A.15) 



The second term in the parentheses appears when the splitting flA.lSP is applied for the 
calculation of 

Using an expansion similar to ( lA.Sp 

F(e"^) — )■ - cos7ra[ln20 — i— 1 sinvro; 

and rescaling the variable (p = t6, one concludes that 



(A.16) 



V{6) 



Re 



t^ 



'Kt{t - i 

The remaining integrals are calculated as above: 



1 TT 

- COS 7rQ;[ln(2t5) — i— ] 



7r 



sm TTO; 



dt 



(A.17) 



t" 



-dt 



177 



-i7ra/2 



t(t — i) sm vra 
Using these values, one gets 



t" Int 



-dt 



ITT 



sm vra 



'ina/2 



cot TTO; 



As T(l - 5) — ^ (2/5)" we conclude that 

(S-i-O 



-- — cot ira ) .(A.l^ 



(A.19) 



TT 



'^(1-'^)= 7 [e(^(l)- J(0)) + l] + -cot(7ra)e + O(5^-") . (A.20) 



4 
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